Quantum effects on the dynamics of a two-mode atom-molecule Bose-Einstein 

condensate 



A. Vardi 1 , V. A. Yurovsky 1 ' 2 , and J. R. Anglin 3 

1 ITAMP, Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138 
2 School of Chemistry, Tel Aviv University, 69978 Tel Aviv, Israel 
3 Center for Ultracold Atoms, MIT 26-237, 77 Massachusetts Avenue, Cambridge MA 02139 



o : 
o . 

(N ■ 

. 



(N 
(N 



> 
On 
CO 

in 
o 



i 

C 

o 
o 



We study the system of coupled atomic and molecular con- 
densates within the two-mode model and beyond mean-field 
theory (MFT). Large amplitude atom-molecule coherent os- 
cillations are shown to be damped by the rapid growth of 
fluctuations near the dynamically unstable molecular mode. 
This result contradicts earlier predictions about the recovery 
of atom-molecule oscillations in the two-mode limit. The fre- 
quency of the damped oscillation is also shown to scale as 
\/~N / log N with the total number of atoms N, rather than 
the expected pure y/N scaling. Using a linearized model, we 
obtain analytical expressions for the initial depletion of the 
molecular condensate in the vicinity of the instability, and 
show that the important effect neglected by mean field theory 
is an initially non-exponential 'spontaneous' dissociation into 
the atomic vacuum. Starting with a small population in the 
atomic mode, the initial dissociation rate is sensitive to the 
exact atomic amplitudes, with the fastest (super-exponential) 
rate observed for the entangled state, formed by spontaneous 
dissociation. 



Recent photoassociation jjj and Feshbach resonance 
[^|,^) experiments suggest the possibility of producing 
molecular Bose-Einstein condensates (BEC) Large 
amplitude coherent oscillations between an atomic BEC 
and a molecular BEC have been theoretically predicted 
||-|| . A common theme to these studies is the use of the 
Gross-Pitaevskii (GP) mean-field theory (MFT), reduc- 
ing the full multi-body problem into a set of two coupled 
nonlinear Schrodinger equations. These are then solved 
numerically to obtain the Josephson-type dynamics of 
the coupled atomic and molecular fields. 

The simple GP dynamics is substantially affected 
by condensate depiction due to inelastic collisions 
P,[7|JlO|l, spontaneous emission, and the inclusion of non- 
condensate modes pO|-p^|. Consequently, the atom- 
molecule oscillations are expected to be strongly damped 
under current experimental conditions. The proposed 
remedy for this detrimental effect involves a recovery of 
an effective two- mode dynamics [ 113] , thereby preventing 
the buildup of thermal population. 

In this article we point out that even in the perfect two- 
mode limit, MFT fails to provide long-term predictions 
due to strong interparticlc entanglement near the dynam- 
ically unstable molecular mode. Quantum corrections to 
MFT appear in the vicinity of its dynamical instabili- 
ties, on time scales that grow only logarithmically with 
the number N of condensate particles [I15H17J- Thus, 



even in the absence of a 'proper' thermal bath of non- 
condensate states, the mean-field equations are coupled 
to a reservoir of Bogoliubov fluctuations pq , [18| . The 
rapid growth of these fluctuations near the instability is 
analogous to the rapid population of the thermal cloud, 
similarly inhibiting the mean-field atom-molecule oscil- 
lations. Our results, obtained using the numerical solu- 
tion of exact quantum equations, go beyond the Hartree- 
Fock-Bogoliubov approach [Q . The leading quantum ef- 
fect is identified as a non-exponential spontaneous decay 
of the molecular condensate and the dynamics is shown 
to be highly sensitive to the initial conditions. We note 
that similar quantum corrections have been predicted for 
parametric oscillations in quantum optics |p.9| . 

We consider the simplest model of the atom-molecule 
condensate, in which particles can only populate two 
second-quantized modes: an atomic mode, associated 
with the creation and annihilation operators a and a+ 
and a molecular mode, associated with the creation and 
annihilation operators b and 6+ . The two modes are cou- 
pled by means of a near-resonant two-photon transition 
or a Feshbach resonance, with a coupling frequency fi 
and detuning A. Setting the zero energy to the energy 
of the molecular mode the two-mode Hamiltonian reads 
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For A = 0, the Hamiltonian of Eq. (|T]) is identical [Q] to 
the well known Hamiltonian describing the optical pro- 
cess of parametric oscillations j20(| , where dissociation is 
equivalent to parametric downconversion and association 
is the analog of second-harmonic generation. We will take 
fl to be real and positive without loss of generality, since 
the relative phase between the modes is determined up 
to an additive constant and the overall sign of H is in- 
significant. 

We obtain a generalization of the Bloch representation 
for the two-mode system (similarly to the approach taken 
in Refs. (l^Jl^]), by introducing the three operators, 



t x = V2 



a'a'b + b'aa 

7V3/2 

a^a^b — b^aa 



(2) 



C z = 



26+6 — fita 
N 



where N denotes the total number of atoms. Discarding 
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c-number terms, the Hamiltonian of Eq. ([[]) then takes 
the simple form 



ihp = [H, p] , 



(J) 



H = h 



N\ 1 * N ~ 

— nc x - —ac z 
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(3) 



and the Heisenberg equations of motion for the three op- 
erators of Eq. (0) read 



dt x 
dt Cy 



-AC,, 



dt 



C z = 




-v ■ 



These three operators do not represent SU(2); but all 
three commute with the conserved total atom number 
N = a) a + 2&t£. 

The mean field approximation is invoked by approxi- 
mating second order expectation values (CiCj) as prod- 
ucts of the first order moments (£j) and (Cj) p!q , ^7| 



(A4->« (A)(4) • 



(•5) 



Approximation (0) is correct to 0(l/y/~N). Thus we can 
also neglect the c-number term y/2/Nil in Eq. (Q). 
Defining s = ((£ x ) , (£ y ) , (£ z )), rescaling the time as 
r = \Z~NClt and using Eq. (g), we obtain the mean-field 
equations 



dT Sy 



dr 
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(6) 



where the dimensionless rescaled detuning 8 is defined 
as 5 = A/y/Nfi, (Eqs. (§) are equivalent to Eqs. (32) of 
Rcf. without the inelastic collision terms). Lyapunov 
analysis of Eqs. @ shows that as long as \6\ < y/2 the 
stationary point s = (0,0, 1) corresponding to the entire 
population being in the molecular mode, is dynamically 
unstable because any small perturbation to the mean- 
field equations (0) would trigger the parametric oscilla- 
tion. In the vicinity of this point, MFT is expected to 
break down on a time scale which is only logarithmic in 
N [lTj . In order to verify this prediction, we solve 
the full N-body problem by fixing N, thereby restrict- 
ing the available phase-space to Fock states of the type 
|n, (N — n)/2) with n atoms and (N — n)/2 molecules, 
where n = 0, 2, 4, N. Thus we obtain an N/2 + 1 
dimensional representation for the Hamiltonian and the 
iV-body density operator p. The quantum solution is 
then obtained by numerically propagating p according to 
the Liouville von-Ncumann equation 



and the expectation values of the three operators of Eq. 
(||) are retrieved as Sj = Tr(pCi). 

In what follows, we shall assume that 5 — 0, as required 
to obtain unity amplitude atom- molecule oscillations. In 
Fig. 1 we plot the expectation value s z , corresponding to 
the population difference between the modes, as a func- 
tion of the rescaled time r, for various values of the total 
particle number. The initial conditions are s = (0, 0,-1), 
corresponding to an initially populated atomic mode. 
The mean field solution (dotted line) depicts the con- 
vergence of s z to the unstable fixed point. The quan- 
tum solutions (identical to the results of Ref. Q obtained 
by solving the iV-particle Schrodinger equation) , initially 
follow the mean-field evolution closely. However, in the 
vicinity of the molecular mode (s z = 1) the quantum tra- 
jectories break away from the mean-field prediction on a 
timescale that grows only logarithmically with N. Thus, 
the oscillation frequency scales with v^V / log(iV) as op- 
posed to the expected scaling with y/~N. Moreover, the 
oscillations are damped by the strong entanglement near 
the molecular mode, in complete analogy with the damp- 
ing of oscillations when the two-mode system is coupled 
to external thermal modes. Full Rabi-type coherent os- 
cillations can only be observed for a single pair of atoms. 

We note that the results of Fig. 1 can not be re- 
produced by a Hartree-Fock-Bogoliubov approach jfj). 
In order to obtain the damping of coherent oscillations, 
one has to go deeper in the Bogoliubov-Born-Green- 
Kirkwood-Yvon (BBGKY) hierarchy of equations of mo- 
tion, and maintain a number of equations comparable to 
the total number of particles N. 
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FIG. 1. Numerically calculated population imbalance be- 
tween atoms and molecules, as a function of the rescaled time 

r for 2 ( ), 10 ( ), 100 (- - -), and 1000 ( ) particles, 

compared with the mean-field prediction (•■•)• 
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Equations (||) are equivalent to the equations of motion 
for the two annihilation operators a and b, 



■ A „ 
ia= —a 4 

ib— — aa 
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Qba f 



(8a) 
(8b) 



To appreciate why MFT fails as it does, we will now focus 
our attention on the vicinity of the dynamically unsta- 
ble all-molecule state, using a linearized model in which 
the molecular annihilation operator b is replaced by a c- 
number 6 of 0(y/N/2). In this approximation, which is 
valid as long as the population of the molecular state is 
large and the effect of its depletion on the atomic popu- 
lation growth rate can be neglected, Eq. (^) becomes 



i a 



(9) 



Equation (g) in combination with its complex conjugate 
provide an autonomic set of two linear operator equa- 
tions, which can be solved using common methods re- 
ducing it to an eigenproblem. When the molecular mode 
is dynamically unstable (\6\ < \/2) we have 



Q\b\ > |A|/2 , 



(10) 



and the exact solution of Eq. takes the rapidly grow- 
ing form 



a(t) = a(0) cosh Xt 



-a{0) + nba\0) 



sinh Xt 



where 



A = ^ 2 H 2 ~(A/2) 2 



(11) 



(12) 



The time-dependence of the atom number operator 
n a (t) = a'(i)a(t) is thus given as 



(n a (t)} = n sp (t) + n st (t) , 
where the term 



n sp {t) 



— ^ sinh 2 Xt 
X 1 



(13) 



(14) 



not accounted for by MFT, corresponds to spontaneous 
dissociation into the atomic vacuum and the term 



n st (t) = (ft o (0)) ( cosh2Ai -I- ^ sinh 2 Xt 



n 
x 



-Im 



2A 2 

j sinh2A/ - iy sinh 2 Xt 



(15) 



depicts stimulated dissociation taking place when the 
atomic state is initially populated. The two terms on 



the r.h.s of Eq. ( p_5[ ) correspond to non-coherent and co- 
herent initial occupation, respectively. 

Starting from the dynamical instability (|6| = y/N/2) 
with zero atomic occupation, and assuming zero detun- 
ing, the initial evolution of the expectation value s z (t) is 
given according to Eqs ([l3|)-([l5j) by the purely sponta- 
neous process, 



S ,W = l-ysinh 2 (r/v / 2) 



(16) 



According to Eq. ( |iq ) the initial decay of the atomic 
mode is non-exponential as the leading decay term is 
quadratic rather than linear in t. This behavior is in ac- 
cordance with the initial non-exponential decay of a gen- 
eral spontaneous emission process 21 . The quadratic 
gain of the atomic population is confirmed by the re- 
sults of Fig. 2 where we compare the initial depletion 
of the molecular mode according to Eq. (|l6| ) with exact 
quantum results obtained for various values of N. The 
agreement is initially excellent until the occupation of the 
atomic mode becomes significant compared to the molec- 
ular occupation. Moreover, it is evident from Eq. ( |l6| ) 
and confirmed by the results of Fig. 2, that the time at 
which the quantum spontaneous emission term will be- 
come significant, grows only logarithmically with N, in 
agreement with our prediction. 
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FIG. 2. Depletion of the purely molecular mode for 

N = 10, 100 and 1000 particles, according to Eq. @ ( ), 

compared with 'exact' numerical results ( ). 

Equation (|l4|) was already obtained in Ref . , as the 
asymptotic expression at t — * oo. However, at this limit, 
the atomic occupation becomes comparable to the molec- 
ular population, and the depletion of the molecular mode 
should be taken into account. Nevertheless, Eq. ( |l4| ) is 
an exact solution to the model of Eq. ([)]) , applicable at 
small t. 
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It is interesting to note that when |<5| > \[2 the molec- 
ular mode is stabilized. Consequently, the exact solution 
of Eq. A) when 



Sl\b\ < |A|/2 
is an oscillatory function of the form 



a(t) =a(0)cos|A|<- — 
|A| 



where now 



^a(O) + Qbat(0) 



|A| = v/(A/2)2-fi2| 6 |2. 



(17) 

sin \ X\t 
(18) 

(19) 



It is also worth noting that as shown in Fig. 1, the evo- 
lution of a single molecule, coupled to a single atomic 
mode, is alwa ys d escribed by an oscillating solution, sim- 
ilar to Eq. ( |l8|) since Bose enhancement, depicted by 
the exponential gain of Eqs. (|ll|) and (|l^) is a collective 
effect, analogous to lasing. 
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FIG. 3. Depletion of the molecular mode for three 
100 particle states with the same initial value of 
«,(0) = 74/75: (a) (2/3)^ 2 \0, N/2) + i(l/3) 1/2 \2, N/2 - 1), 
(b) (2/3) 1/2 |0, 7V/2) -i(l/3) 1/2 |2,iV/2- 1), and (c) the state 
created by the spontaneous process starting from the pure 

molecular mode, according to Eq. (J^) ( ). Corresponding 

solid lines are 'exact' numerical results. 



formed by spontaneous dissociation, always gains faster 
{n sp (t)/n sp (t) — 2AcothAt > 2A), due to its high en- 
tanglement [\(a(t)a(t))\/(h a (t)) = cothAi > 1], In Fig. 
3 we compare the predictions of Eq. (^) with ex- 
act quantum results for three initial amplitudes, corre- 
sponding to the same initial atomic occupation of 4/3 N. 
The (a) and (b) curves correspond to states of the form 
c |0, N/2) +ic 2 \2, N/2-1) with (c = ^/2/3, c 2 = i^/T/S) 
and (co = y/2/3, c 2 — — i^/l/3), respectively, giving 
(b* (a(0)a(0)) = ±i|6|(n a (0)) (in equivalence to coher- 
ent states). The (c) lines depict the continued decay of 
a squeezed state with the same (n a (0)), formed by the 
purely spontaneous process, starting at r = —0.53. As 
expected, an initial exponential atomic loss is observed 
when C2 = iy/1/3 and an exponential atomic gain is ob- 
served when C2 = —iy/l/3. The gain of the sponta- 
neously produced squeezed state, depicted by Eq. (H), 
is initially non-exponential, becoming so only at later 
times. 

In conclusion, due to the dynamical instability of the 
molecular mode, the atomic ensemble produced by the 
dissociation stage of the parametric oscillation, is highly 
entangled. Consequently, atom-molecule coherent oscil- 
lations are predicted to be damped even in the two-mode 
limit. The leading quantum correction is an initially 
non-exponential spontaneous decay, becoming significant 
on a timescale which only grows as log TV. Stimulated 
processes taking place when there is an initial atomic 
population, are sensitive to different initial conditions 
corresponding to the same s z (0). The fastest (super- 
exponential) decay rate is obtained for the initial state 
formed by spontaneous dissociation. This dependence 
originates in the parametric oscillation phase, driving it 
towards the molecular mode or away from it. It has 
significant implications on the case of multiple atomic 
modes, suggesting that the rapid growth observed when 
these modes are initially partially populated Q would 
be sensitive to the exact atomic amplitudes. 
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Finally, we consider the dependence of the stimulated 
decay term as described by Eq. (15), on the exact ini- 
tial conditions. Starting with an initially small coher- 
ent atomic amplitude [| (a(0)a(0)) | = (n a (0))] and as- 
suming exact resonance (A = 0), the initial evolution 
of n st can be varied from exponential gain to exponen- 
tial decay [n st (t) ~ exp(±Ai)] by controlling the rel- 
ative phase between the atomic and molecular modes 
[(a(0)a(0)) and b = (b)}. Moreover, the squeezed state, 
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